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RESUMEN 



ABSTRACT 

Using the abundance matching technique, we infer the local stellar and bary- 
onic mass-halo mass (M s -Mh and Mb-Mh) relations separately for central blue and 
red galaxies (BGs and RGs). The observational inputs are the SDSS central BG 
and RG Stellar Mass Functions and the measured gas mass-M s relations. For halos 
associated to central BGs, the distinct ACDM Halo Mass Function is used and set 
up to exclude: (i) the observed group/cluster mass function and (ii) halos with a 
central major merger at rcsdshifts z < 0.8. For central RGs, the complement of this 
mass function to the total one is used. At Mh > 10 115 M Q , the M s of RGs tend to 
be higher than those of BGs for a given Mh, the difference not being larger than 1.7. 
At M h < 10 n - 5 M Q , this trend is inverted. For BGs (RGs): (a) the maximum value 
of f s = M s /M h is 0.02ltom9 (°- 034 -o.oi5) and [t is attained at log(M h /M )= 12.0 
(=11.9); (b) f s oc Mh {f s oc M^ 3 ) at the low-mass end while at the high-mass end, 
f s cx M h ~ 0A (f s oc M h ~ - 6 ). The baryon mass fractions, f h =M h /M hl of BGs and 



RGs reach maximum values of / b = 0.028±£;^ and / b = 0.034+^ 



, +0.02.5 



014 , respectively. 



At Mh < 10 M Q , the dependence of f^ on Mh is much steeper for RGs than 
for BGs. We discuss on the differences found in the / S -Mh and /b-Mh relations 
between BGs and RGs in the light of semi-empirical galaxy evolution inferences. 
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1. INTRODUCTION 
The galaxy stellar and baryonic mass functions 
(GSMF and GBMF ', respectively), inferred from the 
observed luminosity function and gas fraction-stellar 
mass (/ g — M s ) relation, contain key statistical in- 
formation to understand the physical processes of 
galaxy formation and evolution. Within the context 
of the popular A Cold Dark Matter (ACDM) hierar- 
chical scenario, dark matter halos are the sites where 
galaxies form and evolve (White & Rees 1978; White 
& Frenk 1991). Hence, a connection between GBMF 
or GSMF and the halo mass function (HMF) is ex- 
pected. The result of such a connection is the galaxy 
stellar and baryonic mass-halo mass relations, M s - 
Mh and Mb-Mh, and their intrinsic scatters, both set 
by complex dynamical and astrophysical processes 
intervening in galaxy formation and evolution (see 
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for recent reviews Baugh 2006; Avila-Reese 2007; 
Benson 2010). In this sense, the M b /M h and M s /M h 
ratios quantify the efficiency at which galaxy and 
star formation proceeds within a halo of mass Mh. 
Therefore, the empirical or semi-empirical inference 
of the Mb-Mh and M s -Mh relations and their scat- 
ters (locally and at other epochs) is nowadays a chal- 
lenge of great relevance in astronomy. 

For simplicity, in statistical studies like those re- 
lated to the GSMF, galaxies are labelled by their 
mass alone. However, by their observed proper- 
ties, correlations, and evolution, galaxies show a 
very different nature, at least for the two major 
groups in which they are classified: the rotationally- 
supportcd disk star-forming (late-type) and the 
pressure-supported spheroid quiescent (early- type). 
In the same way, the evolution of galaxies is ex- 
pected to differ if they are centrals or satellites. The 
main intrinsic processes of galaxy evolution are as- 
sociated to central galaxies, while satellite galax- 
ies undergo several extra astrophysical processes be- 
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cause of the influence of the environment of the cen- 
tral galaxy/halo system in which they were accreted. 
Hence, if the M^-M^ or M s -Mh relations are used 
for constraining galaxy formation and evolution pro- 
cesses, these relations arc required by separate for at 
least the two main families of late- and early-type 
galaxies and taking into account whether the galaxy 
is central or satellite. Fortunately, in the last years 
there appeared several studies, in which a decom- 
position of complete GSMFs by color, concentration 
or other easily measurable indicators of the galaxy 
type was carried out (e.g., Bell et al. 2003; Shao 
et al. 2007; Bernardi et al. 2010). Evenmore, in 
a recent work, Yang, Mo & van den Bosch (2009, 
hereafter YMB09) used the Sloan Digital Sky Sur- 
vey (SDSS) data for obtaining the GSMFs of both 
central and central + satellite galaxies separated in 
each case into blue and red objects. 

With the coming of large galaxy surveys, a big 
effort has been done in constraining the z ~ to- 
tal M s -Mh relation (i) directly by estimating halo 
masses with galaxy-galaxy weak lensing, with kine- 
matics of satellite galaxies or with X-ray studies; and 
(ii) indirectly by linking observed statistical galaxy 
properties (e.g., the galaxy stellar mass function 
GSMF, the two-point correlation function, galaxy 
group catalogs) to the theoretical HMF (see for 
recent reviews and more references Moster et al. 
2010; Behroozi, Conroy & Wechsler 2010, hereafter 
BCW10; More et al. 2011). While the latter ap- 
proach does not imply a measure-based determina- 
tion of halo masses, it is simpler from a practical 
point of view, as it allows to cover larger mass ranges, 
and can be extended to higher rcdshifts than the for- 
mer approach (see recent results in Conroy & Wech- 
sler 2009; Moster et al. 2010; Wang & Jing 2010; 
BCW10). Besides, both the weak lensing and satel- 
lite kinematics methods in practice are (still) statisti- 
cal in the sense that one needs to stack large number 
of galaxies in order to get sufficient signal-to-noise. 
This introduces a significant statistical uncertainty 
in the inferred halo masses. 

The indirect approach for linking galaxy and halo 
masses spans a large variety of methods, among them 
the Halo Occupation Distribution (Peacock & Smith 
2000; Bcrlind & Weinberg 2002; Kravtsov et al. 
2004) and the Conditional Luminosity Function for- 
malisms (Yang et al. 2003, 2004). These formalisms 
introduce a priori functional forms with several pa- 
rameters that should be constrained by the observa- 
tions. Therefore, the final inferred M s -Mh relation is 
actually model-dependent and yet sometimes poorly 
constrained due to degeneracies in the large num- 



ber of parameters. A simpler and more empirical 
method -in the sense that it uses only the GSMF 
(or luminosity function) as input and does not re- 
quire to introduce any model- has been found to 
give reasonable results. This indirect method, called 
the abundance matching technique (hereafter AMT; 
e.g., Marinoni & Hudson 2002; Vale & Ostriker 2004; 
Conroy et al. 2006; Shankar et al. 2006; Conroy & 
Wechsler 2009; Baldry, Glazebrook & Driver 2008; 
Guo et al. 2010; Moster et al. 2010; BCW10), is 
based on the assumption of a monotonic correspon- 
dence between M s and M^; in the limit of zero scat- 
ter in the M s -Mh relation, the halo mass Mh cor- 
responding to a galaxy of stellar mass M s , is found 
by matching the observed cumulative GSMF to the 
theoretical cumulative HMF. 

In this paper we apply the AMT in order to in- 
fer the local M s -Mh relation for central blue and red 
galaxies separately, which requires as input both the 
observed central blue and red GSMFs, taken here 
from YMB09. Note that in order to infer the M s -M h 
relation of galaxy subpopulations (e.g., blue/red or 
central/satellite ones) solely from the overall GSMF, 
models for each subpopulation should be introduced, 
which largely increases the uncertainty in the result. 
Regarding the HMFs to be matched with the cor- 
responding observed central GSMFs, the theoretical 
HMF is decomposed into two functions -associated 
to halos hosting blue and red galaxies- based on em- 
pirical facts: blue galaxies are rare as central objects 
in groups/ clusters of galaxies, and they should not 
have undergone late major mergers because of the 
dynamical fragility of disk (blue) galaxies. Nowa- 
days, it is not clear whether the M s -Mh relation 
varies significantly or not with galaxy color or type. 
Previous studies that discussed this question were 
based on direct methods: the weak lensing (Mandel- 
baum et al. 2006) and satellite kinematics (More et 
al. 2011) techniques. The uncertainties in the results 
of these studies are yet large, and can be subject to 
biases intrinsic to the sample selection and to effects 
of environment. 

We also estimate here the galaxy baryon mass- 
halo mass relations, Mb-Mjjj, where M^= M s + M g , 
by using the GSMFs combined with average ob- 
servational determinations of the galaxy gas mass, 
M g , as a function of M s . The galaxy baryonic mass 
fraction, /b= M^/M^, and its dependence on mass 
is important for constraining models and simula- 
tions of galaxy evolution, and is also a key input 
for some approaches, implemented for modelling the 

4 We assume that the galaxy baryonic mass is included in 
the halo (virial) mass M^. 
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most generic population of galaxies, namely isolated 
(central) disk galaxies (e.g., Mo, Mao & White 1998; 
Firmani & Avila-Reese 2000; van den Bosch 2000; 
Stringer & Benson 2007; Dutton et al. 2007; Gnedin 
ct al. 2007; Dutton & van den Bosch 2009). In these 
and other studies, it was shown that several disk 
galaxy properties, correlations and their scatters de- 
pend (or are constrained) by /b- In a similar way, 
the fb-Mh dependence is expected to play some role 
in the results of structural and dynamical models of 
spheroid-dominated galaxies. 

In Section 2 we describe the method and the data 
input. The stellar/baryon mass-halo mass relations 
for the total, blue and red (sub)samples are pre- 
sented in Section 3. In Section 4 we compare our 
results with other observational works, and discuss 
whether they are consistent or not with expectations 
of semi-empirical inferences. The summary and our 
conclusions are given in Section 5. 

2. THE METHOD 
The AM statistical technique is based on the hy- 
pothesis of a one-to-one monotonic increasing rela- 
tionship between M s (or Mb) and Mh. Therefore, 
by matching the cumulative galaxy stellar and halo 
mass functions, for a given M s a unique Mh is as- 
signed: 

/>oo />oo 

/ 4> h {M' h )dM' h = / MK)dM' s , (1) 
JM h JM S 

where 4>h is the overall HMF (distinct + subhalos) 
and c/> s is the overall GSMF; distinct halos are those 
not contained inside more massive halos. It is rea- 
sonable to link central galaxies with distinct halos. 
Therefore, in the case of using the GSMF for only 
central galaxies, the distinct HMF should be used for 
the matching. Since the main purpose of this paper 
is the inference of the M s -Mh (and the correspond- 
ing Mb-Mh) relation for blue (red) galaxies, (i) a 
GSMF that separates galaxies by color is necessary 
(the data to be used here are discussed in §§2.1), and 
(ii) a criterion to select the halos that will likely host 
blue (red) galaxies shall be introduced (see §§2.2.1). 

In this paper we will not carry out an exhaus- 
tive analysis of uncertainties in the inference of the 
M s -Mh relation with the AMT. This was extensively 
done in BCW10 (see also Moster et al. 2009). In 
BCW10 the uncertainty sources are separated into 
three classes: uncertainties (i) in the observational 
inference of GSMF, (ii) in the dark matter HMF, 
which includes uncertainties in the cosmological pa- 
rameters, and (iii) in the matching process arising 
primarily from the intrinsic scatter between M s and 
M h . 



2.1. Galaxy and Baryonic Stellar Mass Functions 

In the last years, complete galaxy luminosity 
functions (and therefore, GSMFs) were determined 
for local samples covering a large range of lumi- 
nosities (masses). The stellar mass is inferred from 
(multi)photometric and/or spectral data (i) by us- 
ing average stellar mass-to-light ratios, depending 
only on color (inferred from application of stellar 
population synthesis -SPS- models to galaxy sam- 
ples with independent mass estimates, e.g. Bell et 
al. 2003), or (ii) by applying directly the SPS tech- 
nique to each sample galaxy, when extensive multi- 
wavelength and/or spectral information is available. 

In both large uncertainty is introduced 

in the inference of M s due to the uncertainties in the 
IMF, stellar evolution, stellar spectral libraries, dust 
extinction, metallicity, etc. Bell et al. (2003) esti- 
mated a scatter of « 0.1 dex in their M s /L ratios 
in infrared bands. Conroy, Gunn & White (2009) 
carried out a deep analysis of propagation of uncer- 
tainties in SPS modelling and concluded that M s at 
z ~ carry errors up to ~ 0.3 dex (but see Gal- 
lazzi & Bell(2009)). Here, we will consider an over- 
all systematical uncertainty of 0.25 dex in the M s 
determination (see BCW10). 

Most of the current local GSMFs were inferred 
from 2dF Galaxy Redshift Survey, Two Micron All- 
Sky Survey (2MASS) and SDSS (e.g., Cole et al. 
2001; Bell ct al. 2003; Baldry et al. 2006). The 
low-mass completeness limit due to missing of low 
surface brightness galaxies is at ~ 10 8 5 M s (Baldry 
ct al. 2008). An upturn of the GSMF close to this 
end (below M s ^ 10 9 M Q ) was confirmed in several 
recent works (Baldry et al. 2008; YMB09; Li & 
White 2009). Due to this upturn, a better fit to the 
GSMFs is obtained by using a double or even triple 
Schechter function. Since the low-mass end of the 
GSMF is dominated by late- type galaxies, this up- 
turn plays an important role in the M s -Mh relation 
of late- type galaxies at low masses. 

For our purposes, observational works where the 
GSMF is decomposed into late- and early-types 
galaxies are required. Such a decomposition has been 
done, for example, in Bell et al (2003), who combined 
22679 SDSS Early Data Release and 2MASS galax- 
ies, and used two different criteria, color and concen- 
tration, to split the sample into two types of galaxies. 
A much larger sample taken from the NYU-VAGC 
based on the SDSS DR4 has been used by YMB09 
(see also Yang, Mo & van den Bosch 2008), who split 
the sample into blue and red sub-samples according 
to a criterion in the °- 1 (g — r) — M r diagram. In both 
works, M s is calculated from the r— band magnitude 
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Fig. 1. Left panel: Different local GSMFs for all galaxies. The reported data in Bell et al. (2003, pink squares) and 
YMB09 (black hexagons) are plotted directly, while for Baldry et al. (2008, blue solid line) and Li & White (2009, dot- 
dashed green line), the best fits these authors find to their samples are plotted. Red triangles show the data from YMB09 
corresponding to the GSMF of central-only galaxies. Middle and left panels: Data corresponding to the decomposition 
of the GSMF into blue and red galaxies, respectively, from Bell et al. (2003) and for the all and central-only galaxies 
from YMB09. 



by using the corresponding color-dependent M s /L r 
ratio given in Bell et al. (2003). In YMB09 each 
color subsample is in turn separated into central and 
satellite galaxies according to their memberships in 
the constructed groups, where the central galaxy is 
defined as the most massive one in the group and the 
rest as satellite galaxies. 

In Figure [Q the Bell et al. (2003) and YMB09 
GSMFs arc reproduced by using the data sets re- 
ported in these papers. In the left panel, the full 
sample from each work (solid squares and solid 
hexagons, respectively) are plotted, as well as the 
case of central-only galaxies from YMB09 (solid tri- 
angles); both GSMFs and the other ones plotted 
in this figure are normalised to h = 0.7 and to a 
Chabricr (2003) IMF. In the central and right panels, 
the corresponding blue (late-type) and red (early- 
type) sub-samples are plotted with the same symbols 
on the left panel. For the Bell et al. sub-samples, 
only those separated by their color criterion arc plot- 
ted. Both GSMFs corresponding to the full and blue 
sub-samples are in good agreement for M s ^ 10 95 
Mq. For lower masses, the Bell et al. GSMF's are 
higher. On one hand, the Bell et al. sample is much 
smaller than the YMB09 one (therefore its cosmic 
variance is more significant). On the other hand, 



the redshift completeness and M s limit in YMB09 is 
treated with updated criteria. 

In Figure [TJ we also plot fits to the overall GSMF 
presented in Baldry et al. (2008, double Schechter 
function, solid blue line) and in Li & White (2009, 
triple Schechter function, dashed green line) for new 
SDSS releases and by using directly SPS models 
to estimate M s for each galaxy. These fits agree 
well with the YMB09 data in the mass range 9.2 ^ 
log(Af s /M Q ) ^ 11.2. For smaller masses, the Baldry 
et al. fit tends to be steeper while the Li & White 
fit tends to be shallower than the YMB09 data. For 
larger masses, both fits decrease faster with M s than 
the YMB09 data. All these (small) differences are 
due to the different methods used to estimate M s , as 
well as the different volumes and limit corrections of 
the samples (sec Baldry et al. 2008, YMB09, and Li 
& White 2009 for discussions). 

The split into two colors of the sample used by 
YMB09 is a rough approximation to the two main 
families of disk- and spheroid-dominated galaxies. 
It is well known that the morphological type corre- 
lates with the galaxy color, though with a large scat- 
ter. There is for example a non- negligible fraction of 
galaxies (mostly highly inclined) that are red but of 
disk-like type (e.g., Bernardi et al. 2010). However, 
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Fig. 2. Gas mass vs stellar mass for a sample of disk 
high and low surface brightness galaxies collected and 
homogenised by Avila-Reese et al. (2008, blue dots with 
error bars) and for a sample of disk galaxies presented by 
McGaugh (2005, blue crosses). The solid blue line is the 
orthogonal linear doubly-weighted regression to the data 
from the former authors and the dashed lines show an 
estimate of the intrinsic scatter around the fit. The solid 
red line is an estimate of the M g -M s correlation for red 
galaxies by using our fit to blue galaxies and the ratio of 
blue-to-red atomic gas fraction determined in Wei et al. 
(2010; see text). 

given that here we consider a partition of the overall 
sample just in two groups, we believe that is rea- 
sonable to assume in a first approximation that the 
color criterion for the partition will provide a similar 
result at this level as a morphological criterion. 

For the YMB09 sample, the blue and red galaxies 
are 55% and « 45%, respectively, for M s > 3 x 10 s 
M Q . Red galaxies dominate the total GSMF at large 
masses. At M s « 2 x 10 10 M Q the abundances of red 
and blue galaxies are similar and at lower masses the 
latter are increasingly more abundant than the for- 
mer as M s is smaller. For AI S 10 9 M Q , the abun- 
dance of red galaxies, mainly central ones, steeply 
increases towards smaller masses. The existence of 
this peculiar population of faint central red galaxies 
is discussed in YMB09. Wang et al. (2009) sug- 
gested that these galaxies are hosted by small halos 
that have passed through their massive neighbors, 
and the same environmental effects that cause satel- 
lite galaxies to become red are also responsible for 
the red colors of such galaxies. However, as these au- 
thors showed, even if the environmental effects work, 
there are in any case over 30% of small halos that are 
completely isolated in such a way that these effects 
can not be invoked for them. 

In the YMB09 sample, around 70% of the galax- 



ies are central. As mentioned in the Introduction, 
the inference of the M s -Mh relation for central-only 
galaxies is important for studies aimed to constrain 
galaxy formation and evolution in general; satellite 
galaxies are interesting on its own but they lack gen- 
erality because their evolution and properties are af- 
fected by extra environmental processes. 

In what follows, the YMB09 GSMF provided 
in tabular form and split into blue/red and cen- 
tral/satellite galaxies will be used for applying the 
AMT. Our main goal is to infer the M s -Mh rela- 
tion for central blue (late-type) and red (early-type) 
galaxies. 

We will infer also the corresponding M^-M^ 
(baryonic) relations. The blue and red GBMFs are 
estimated from the blue and red GSMFs, respec- 
tively, where in order to pass from M s to Mb, the 
cool (atomic and molecular) gas mass, M g , corre- 
sponding on average to a given M s is taken from the 
empirical blue and red M g -M s relations. In Figure 
[21 a compilation of observational estimates is plot- 
ted in the M s -M g plane for a sample of disk galaxies 
that includes low surface brightness galaxies from 
Avila-Reese et al. (2008; blue dots with error bars; 
they added H2 mass contribution by using an esti- 
mate for the H2-to-HI mass ratio as a function of 
galaxy type), and for another galaxy sample from 
McGaugh (2005; blue crosses; no H2 contribution is 
considered and their dwarf galaxies were excluded). 
An orthogonal linear doubly-weighted regression to 
the data from Avila-Reese et al. (2008) gives: 

^^- 43 x(^M°' 62 (2) 
10i°M Q -" A6x ^lOiOMj ' {Z) 

This fit is plotted in Figure [2] with its corresponding 
estimated scatter (« 0.3 dex; blue solid and dashed 
lines). This is the relation and its scatter used to 
calculate Mb and the blue GBMF. A similar rela- 
tion has been inferred by Stewart et al. (2009). The 
gas fractions in red galaxies are much smaller than 
in the blue galaxies. For sub-samples of blue and red 
galaxies, Wei et al. (2010) reported for each one the 
atomic gas fractions versus M s (molecular gas was 
not included). The ratio of their fits to these data 
as a function of M s is used here to estimate from eq. 
([2l blue galaxies) the corresponding average M g for 
red galaxies as a function of M s . The red solid line 
shows the obtained relationship. As an approxima- 
tion to the scatter (short-dashed lines), the average 
scatter reported for red galaxies in Wei et al. (2010) 
is adopted here. 
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2.2. Halo and sub-halo mass functions 

A great effort has been done in the last decade to 
determine the HMF at z = and at higher redshifts 
in N-body cosmological simulations. A good fit to 
the results, at least for low redshifts, is the universal 
function derived from a Press-Schechter formalism 
(Press & Schcchtcr 1974) generalized to the ellip- 
tical gravitational collapse (Sheth & Tormen 1999 
hereafter S-T). In fact, Tinker et al. (2008) have 
shown that at the level of high precision, the HMF 
may change for different cosmological models and 
halo mass definitions as well as a function of z. For 
our purposes and for the cosmology used here, the S- 
T approximation provides a good description of the 
z = HMF of distinct halos: 



2 puv 
- Ml 



d In a 



d\nM h 



exp 



<j> h {M h )dM h = 

21 

dM h (3) 



v 

Y 



where A = 0.322, q = 0.3, v 2 = a(5 c / D(z)a{M h )) 
with a = 0.707, 5 C = 1.686fi^ 0055 is the linear 
threshold in the case for spherical collapse in a flat 
universe with cosmological constant, D{z) is the 
growth factor and cr(Mh) is the mass power spec- 
trum variance of fluctuations linearly extrapolated 
to z = 0. The halo (virial) mass, Mh is defined 
in this paper as the mass enclosed within the ra- 
dius where the overdensity is p vlI — A times the 
mean matter density, p~M\ A k, 340 according to 
the spherical collapse model for the cosmology used 
here. The cosmological parameters assumed here are 
close to those of WMAP5 (Komatsu et al. 2009): 
fl M = 0.27, fi A = 1 - Hm = 0.73, h = 0.70, cr 8 = 0.8. 

The distinct HMF should be corrected when a 
GSMF corresponding to all galaxies is used in the 
AMT. In this case, satellite galaxies are included in 
the GSMF. Therefore, subhalos should be taken into 
account in the HMF. The subhalo fraction is not 
more than « 20% of all the halos at z = (e.g., 
Shankar et al. 2006; Conroy et al. 2006; Giocoli et 
al. 2010; BCW10). When necessary, we correct the 
S-T HMF for (present-day) subhalo population by 
using the fitting formula to numerical results given 
in Giaccoli et al. (2010): 



m sub ) 



din m s 



^ m L ! ex P 



m 



(4) 



with i] = 0.07930, logA = 7.812, log(m o /M ) = 
13.10 and 7 = 0.407. 

The upper panel of Figure [3] shows the (distinct) 
S-T HMF (solid line), the sub-halo HMF (short- 
long-dashed line), and the distinct+subhalo HMF 



(dash-dotted line). The correction by sub- halos in 
the abundance is small at low masses and negligible 
at high masses. When the GSMF refers only to cen- 
tral galaxies -which is the case in this paper-, then 
it is adequate to use namely the distinct HMF for 
the AMT, i.e. the subhalo abundance correction is 
not necessary. 

2.2.1. Haloes hosting blue and red galaxies 

In the AMT, the cumulative GSMF and HMF 
are matched in order to link a given M s to M^. 
When a subsample of the total GSMF is used -as 
is the case for inferring the M s -Mh relation of only 
late- or early-type galaxies- it would not be correct 
to use the total HMF for the matching. This func- 
tion, in the ignorance of which is the mass function 
of halos hosting blue (red) galaxies, at least should 
be re-normalised (decreased uniformly) by the same 
fraction corresponding to the decrease of the sub- 
sample GSMF with respect to the total GSMF. In 
YMB09, « 55% (« 45%) of the galaxies are in the 
blue (red) sub-samples for M s > 3 x 10 8 M©. We 
may go one step further by proposing general ob- 
servational/physical conditions for halos to be the 
host of blue (late-type) or red (early-type) galaxies. 
Note that the division we do here of galaxies is quite 
broad -just in two groups-, therefore very general 
conditions are enough. 

Haloes that host central blue and red galaxies 
are expected to have (i) a different environment, and 
(ii) a different merger history. We take into account 
these two factors in order to roughly estimate the 
HMF of those halos that will host today central blue 
and red galaxies. 

Environment.- Blue (late- type) galaxies are rare 
in the centers of groups and clusters of galax- 
ies (high-density environments; e.g., Norberg et al. 
2001; Zehavi et al. 2005; Li et al. 2006; de Lap- 
parent & Slezak 2007; Padilla, Lambas & Gonzalez 
2010; Blanton & Moustakas 2009, and more refer- 
ences therein). For example, in the SDSS YMB09 
sample that we use here (see also Wcinmann et al. 
2006), among the groups with 3 or more members, 
the fraction of those with a central blue galaxy is only 
~ 20%, and most of these central galaxies are actu- 
ally of low masses. Therefore, cluster- and group- 
sized halos (more massive than a given mass) can 
not be associated to central blue galaxies when using 
the AMT. This means that the halo mass function of 
groups/clusters of galaxies should be excluded from 
the theoretical HMF (Shankar et al. 2006). 

Heinamaki et al. (2003) have determined the 
HMF of groups with 3 or more members and with 
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Fig. 3. Upper panel: Distinct S-T HMF for the cos- 
mology adopted in this paper (solid black line) , sub-halo 
mass function at z = according to Giaccoli et al. (2010, 
short- long-dashed purple line), and the sum of both (dot- 
dashed orange line). The solid dots are measures of the 
group/cluster mass function according to Heinamaki et 
al. (2003) and adequately corrected to our definition of 
virial halo mass; the dot-long-dashed cyan line is a eye-fit 
to the data. Lower panel: The same distinct S-T HMF 
(solid black line) shown in the upper panel but (i) ex- 
cluding the halos that suffered late major mergers -since 
z = 0.8- (short-dashed black line) and (ii) excluding 
these halos and those of observed groups/clusters (long- 
dashed blue line). The latter is the HMF to be assigned 
to the sub-sample of central blue galaxies. The comple- 
ment of this function to the total (S-T) one (dot-dashed 
red line) is the HMF to be assigned to the sub-sample of 
central red galaxies. The inset shows the ratio of number 
densities of halos that did not suffer major mergers since 
z = 0.8 to all the (distinct) halos according to measures 
in a cosmological N-body simulation (Colin et al. 2011, 
see text). The fit to this ratio (solid line in the inset) is 
what has been used to correct the S-T HMF for halos 
thad did not suffer late major mergers. 



a number density enhancement Sn/n > 80 from the 
Las Campanas Redshift Survey. The authors esti- 
mated the corresponding group virial mass on the 
basis of the line-of-sight velocity and harmonic ra- 
dius of the group, in such a way that this mass is 
defined at the radius where 5n/n=80. The observa- 
tional galaxy overdensity Sn/n is related to the mass 
overdensity Sp/ p roughly through the bias parameter 
b: Sp/p =(l/b)xSN/N, where b m l/cr 8 (Martinez 
et al. 2002). Hence, for ct 8 = 0.8, Sp/p 64; 
since the group selection was carried out in Tucker 
et al. (2000), where an Einstein-de Sitter cosmology 
was used, then p = p crit in this case. In our case, 
the halo virial mass is defined at the radius where 
Sp/p « 340 (see §§2.2); in terms of p c riti our over- 
density is 340 x fijvf = 92. Therefore, the halo virial 
masses in Heinamaki ct al. (2003) should be slightly 
larger than those used here. For the NFW halos of 
masses larger than ~ lfj 13 M Q , the differences are 
estimated to be factors 1.10-1.20. We correct the 
group masses of Heinamaki et al. (2003) by 15%. In 
the upper panel of Fig. [3l the corrected group (halo) 
mass function is reproduced (solid dots) and a eye-fit 
to them is plotted (dot-dashed cyan line). 

Merger history.- Disk (blue, late-type) galaxies 
are dynamically fragile systems and thus they are 
not expected to survive strong perturbations such as 
those produced in major mergers or close interac- 
tions. However, as several theoretical studies have 
shown (e.g., Robertson et al. 2004; Governato et al. 
2008), when the mergers are gas-rich ('wet') and/or 
at early epochs (in fact, both facts are expected to 
be correlated), it is highly probable that a gaseous 
disk is regenerated or formed again with the late ac- 
creted gas. Therefore, a reasonable restriction for 
halos that will host disk galaxies is that they did not 
undergo central major mergers since a given epoch 
(at earlier epochs, while the central major merger 
may destroy the disk, a new gaseous disk can be 
formed later on). Based on numerical simulations, 
Governato et al. (2008) suggested that a 'wet' ma- 
jor merger of disk galaxies at z ~ 0.8 has yet a non- 
negligible probability of rebuilding a significant disk 
by z ~ 0. We will assume here that halos for which 
their centers have a major merger at z < 0.8 will not 
host a disk galaxy. 

In Colm et al. (2011; in prep.) the present-day 
abundance fraction of halos with no central major 
merger since z = 0.8 was measured as a function 
of Mh from an N-body ACDM cosmological high 
resolution simulation with f2 m = 0.24, Q\ = 0.76, 
and as = 0.75 (box size and mass per particle of 
64 ft, _1 Mpc and 1.64 x 10 7 /i _1 M o , respectively). 
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The friends-of-friends (FOF) method with a linking- 
lcngth parameter of 0.17 was applied for identifying 
halos. The mass ratio to define a major merger was 
q = Afh,2/A^h,i > 0.2 and the merger epoch was 
estimated as that one when the center of the ac- 
creted halo arrived to the center of the larger halo 
by dynamical friction; this epoch is calculated as the 
cosmic time when both FOF halos have "touched" 
plus the respective dynamical friction (merging) time 
as given by the approximation of Boylan-Kolchin et 
al. (2008). The fraction of halos that did not suf- 
fer a major merger since z = 0.8 with respect to all 
the halos as a function of Mh measured in Colin et 
al. (2011) is used here to correct our distinct S-T 
HMF. This measured fraction is showed in the inset 
in the lower panel of Fig. [3J the solid line is a lin- 
ear fit by eye in the log-log plot: log(n no MM/n a u)= 
0.472 - 0.0651og(M h /M Q ). As it is seen, the fraction 
slightly decreases with mass, which is consistent with 
the idea that larger mass halos are assembling later 
with a significant fraction of their masses being ac- 
quired in late major mergers. After the correction 
mentioned above, we get the mass function of halos 
that did not suffer a central major merger (q > 0.2) 
since z = 0.8 (short-dashed black line in the lower 
panel of Fig. [3]). 

The final corrected HMFs.- The function ob- 
tained after (i) subtracting from the distinct S-T 
HMF the group mass function and (ii) excluding ha- 
los that did not suffer a late central major merger is 
plotted in Fig. |3j(bluc long-dashed line). This mass 
function is proposed here to correspond to halos that 
host today blue galaxies. The overall number frac- 
tion of these halos with respect to the distinct ones 
(described by the S-T HMF) is ~ 58%, which is 
roughly consistent with the fraction of blue galax- 
ies in the YMB09 sample. The HMF corresponding 
to the complement is plotted in Fig. [3] as the red 
dot-dashed curve. By exclusion, this HMF will be 
associated with the GSMF of the red central galaxy 
sub-sample for deriving the M s -Mh relation of red 
galaxies. 

3. RESULTS 

3.1. The overall, central, and satellite stellar-halo 
mass relations 

In Fig. |4l the M s -Mh relation obtained by us- 
ing the Li & White (2009) GSMF (see §§2.1 and 
Fig. [J) and the S-T HMF corrected to include sub- 
halos is plotted (long-dashed blue line). The rela- 
tion given by BCW10, who also used as input the 
Li & White (2009) GSMF, is shown (short-dashed 
red line). Both curves are almost indistinguishable, 



showing an excellent consistency of our results with 
those of BCW10 in spite of the differences in some 
of the methodological aspects. 

Further, we plot in Fig. [4] the M s -Mh relation as 
above but using now the total YMB09 GSMF (dot- 
dashed pink line). This relation is similar to the 
one inferred using the Li & White (2009) GSMF. 
For log(Mh/M Q ) ^ 12, the former shifts with mass 
slightly to higher values of M s for a given Mh than 
the latter (at log(M h /M Q ) = 13.5 the difference is 
not larger than 0.08 dex in logAf s ). Such a shift is 
explained by the (small) systematical difference be- 
tween the YMB09 and Li & White GSMFs at masses 
larger than log(M s /M ) - 11 (see §§2.1 and Fig. p. 

In Fig. |U the M s -Mh relations given in Baldry 
et al. (2008, dot-dashed orange line), Moster et al. 
(2010, short-long-dashed line) and Guo et al. (2010, 
dotted green line) are also plotted. When it was nec- 
essary, we have corrected the stellar masses to the 
Chabrier IMF, and the halo masses to the definition 
of virial mass used here (see §§2.2). As mentioned 
above, Baldry et al. corrected their HMF to ex- 
clude groups/clusters of galaxies (something that we 
do but only for the central blue galaxies, see §§2.2.1 
and the result below). As seen in Fig. HI their correc- 
tion produces a steeper M s -Mh relation at the high- 
mass side than in our case. Moster et al. and Guo 
et al. constrained the M s -Mh relation by assigning 
stellar masses to the halos and subhalos of an TV- 
body cosmological simulation in such a way that the 
total GSMF is reproduced. Therefore, by construc- 
tion, their M s -Mh relations take into account the 
group/ cluster halo masses issue. The M s -Mh rela- 
tions in both works are also slightly steeper than ours 
at high masses but shallower on average than that 
one of Baldry et al. (2008). Note that in BCW10 
the scatter in M s at fixed Mh was taken into account 
but the group/cluster halo masses issue was not. 

The M s -M h relation using the YMB09 GSMF 
only for central galaxies and the distinct (S-T) HMF 
is plotted in the lower panel of Fig. 2] (solid black 
line). At large masses, this relation is quite similar 
to the one for all galaxies/satellites and halos/sub- 
halos (dot-dashed pink line) . This is because at large 
masses the great majority of galaxies are centrals and 
the correction by sub-halos is negligible (see Figs. [T] 
and|3j. At lower masses, the exclusion of satellites 
and sub-halos implies a lower M s for a given Mh. 
This is because the GSMF decreases more than the 
HMF as the mass is smaller when passing from the 
total (galaxy and halo) samples to the central-only 
galaxy /distinct halo samples. The physical interpre- 
tation of this result could be that satellite galaxies 
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of a given M s have less massive halos than central 
galaxies due to tidal stripping. The M s -Mh relation 
derived only for the satellites YMB09 GSMF and the 
Giocoli et al. (2010) z = sub-halo HMF is plotted 
in the lower panel of Fig. 2] (short-long-dashed cyan 
line) . 

3.1.1. Uncertainties 

The uncertainty (standard deviation) in the Ms- 
Mh relation obtained by using the YMB09 central 
GSMF and the distinct S-T HMF (solid line), is 
plotted in Fig. @] (grey curves connected by vertical 
lines). As remarked in §2, we did not take into ac- 
count all possible uncertainty sources in the M s -Mh 
relation but have just considered the two following 
ones: 

(i) The systematic uncertainty in stellar mass 
estimates, which is an uncertainty in the GSMF. 
We assume for this uncertainty a scatter of 0.25 
dex (Gaussian distributed) independent of mass, and 
propagate it to the M s -Mh relation (by far it results 
the dominant source of error in the relation obtained 
with the AMT, see below and BCW10). 

(ii) The intrinsic scatter in stellar mass at a fixed 
halo mass, which is an uncertainty in the process 
of matching abundances. To take into account this 
scatter in M s at fixed Mb a probability density distri- 
bution should be assumed. The convolution of this 
distribution with the true or intrinsic GSMF gives 
the measured GSMF. The cumulative true GSMF 
is then the one used for the AM (BCW10). The 
observational data allow to estimate the scatter in 
luminosity (or M s ) and to date it appears to be in- 
dependent of M h (More et al. 2009; Y2009). 
In BCW10 a log-normal mass-independent scatter 
in M s of 0.16±0.04 is assumed. Here, we follow the 
overall procedure of BCW10 for taking into account 
this scatter. 

We also explored the effect of (Hi) the statisti- 
cal uncertainty in the number density of the GSMF 
(as given in YMB09), but we have found that the ef- 
fect is negligible as compared to the one produced by 
item (i) (see also BCW10, their §§4.3.1). The effect 
of the intrinsic scatter in M s for a given Mh is also 
very small in the overall scatter of the M s -Mh rela- 
tion but it affects the high mass end of the calculated 
M s -M h relation, where both the GSMF and HMF 
decay exponentially, since there are more low mass 
galaxies that arc scattered upward than high mass 
galaxies that are scattered downward (BCW10). For 
instance, at Mh = 10 13 5 M Q , the stellar mass af- 
ter including this scatter is 1.2 times smaller. The 
contribution from all other sources of error, includ- 
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Fig. 4. Upper panel: Stellar mass vs halo mass as inferred 
here by using the Li & White (2009) overall GSMF and 
the S-T HMF increased by the subhalo population (long- 
dashed blue line) to be compared with the BCW10 infer- 
ence, who used the same GSMF (short-dashed red line). 
The dot-dashed pink line shows the same M s vs Mh infer- 
ence but by using the overall YMB09 GSMF. Different 
determinations of the overall M s -Mh relation by other 
authors (indicated in the panel), who took into account 
in different ways the issue of group/cluster masses (see 
text) are also plotted Lower panel: Same M a -Mh rela- 
tion as in the upper panel (dot-dashed pink line) but for 
the central-only YMB09 GSMF and the S-T (distinct) 
HMF (solid line). The grey curves connected by verti- 
cal lines show the estimated la uncertainty for the latter 
case. The M s -Mh relation inferred for the only satellite 
YMB09 GSMF and the Giocoli et al. (2010) 2 = sub- 
halo mass function is plotted with the short-long-dashed 
cyan line. 
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ing uncertainties in the cosmological model, is much 
smaller ranging from 0.02 to 0.12 dex at z = 0. 

From Fig. [l]we see that the la uncertainty in the 
M s -Mh relation is approximately 0.25 dex in logM s 
without any systematic dependence on Mh, in good 
agreement with previous results (BCW10; Moster et 
al. 2010). This uncertainty is larger than the differ- 
ences between the M s -Mh average relations found 
by different authors, including those that use the 
indirect AMT but with different GSMFs, method- 
ologies, and corrections, and those who use more 
sophisticated formalisms (see for comparisons and 
discussions BCW10 and More et al. 2011). On one 
hand, this shows that most methods and recent stud- 
ies aimed at relating halo masses to observed galaxies 
as a function of their stellar masses are converging 
to a relatively robust determination. On the other 
hand, this result suggests that attaining a higher pre- 
cision in estimating il/ s from observations is the cru- 
cial task for lowering the uncertainty in the inference 
of the M s -Mh relation. 

3.2. The stellar-halo mass relations for central blue 
and red galaxies 

The upper and lower left panels of Figure [5] show 
the mean M s -M^ and / S -Mh relations for: all cen- 
tral galaxies (solid line, as in Fig. 0]), central blue 
(short-dashed line), and central red (long dashed 
line) galaxies. In order to infer these relations for 
blue galaxies, the central blue YMB09 GSMF and 
the distinct (S-T) HMF corrected for excluding halos 
(i) associated to observed groups/clusters of galax- 
ies and (ii) that suffered central major mergers since 
z = 0.8 (see §§2.2.2) were used. In the case of 
red galaxies, the central red YMB09 GSMF and the 
HMF complementary to the one associated to blue 
galaxies were used. 

The shaded area in Figure [5] is the same la un- 
certainty showed in Fig. @] for the overall central 
sample. The uncertainties corresponding to the M s - 
Mh and / S -Mh relations for the blue and red galaxy 
sub-samples would be close to the one of the total 
sample in case the corrections made to the HMF 
do not introduce an extra uncertainty. In fact this 
is not true, in particular for the group/cluster mass 
function introduced to correct the HMF associated 
to blue galaxies. Unfortunately, the work used for 
this correction does not report uncertainties. Hence, 
the uncertainties calculated here for the blue and red 
samples (shown explicitly in Fig. [5] below) could be 
underestimated, specially at large masses. 

In the mass range 11.5 < log(M h /M Q ) < 13.0, 
the M s -Mh and / S -Mh relations for central blue 
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Fig. 5. Left panels: Mean M a -M h (top) and f s -M h 
(down) relations of all central (solid black line), blue 
central (long-dashed blue line), and red central (short- 
dashed red line) galaxies as inferred here by using the 
YMB09 data. The grey curves connected by vertical lines 
show the lcr uncertainty for the all galaxies case; similar 
uncertainty regions around the main relations are found 
for the blue and red sub-samples (see Fig. [5]). Right 
panels: Same as in left panels but for Mb instead of M s . 
Dotted lines: /b = fu/5 and fu/30, where fu = 0.167 is 
the universal baryon fraction. 
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Fig. 6. Left panels: Mean M s -M h (top) and f s -M h 
(down) relations of central blue (blue lines) and red 
(red lines) galaxies when (i) no systematical corrections 
to the corresponding "blue" and "red" HMFs were ap- 
plied apart from re-normalisations in the global abun- 
dance (see text, solid lines), (ii) the HMFs were cor- 
rected by group/cluster abundances and re-normalised 
(long-dashed lines), and (iii) the HMFs were corrected 
both by group/cluster abundances and late major merg- 
ers (as in Fig. [4] dot-dashed lines). Right panels: Same 
as in left panels but for Mb instead of M s . 



(red) galaxies lie slightly below (above) the rela- 
tions corresponding to the overall sample. For 
masses below these ranges, the trends invert. The 
/ S -Mh curves for blue and red sub-samples peak 
at log(Mh/M Q )= 11.98 and 11.87, with values of 
/„ = 0.021±g:85| and / s = 0.034t°;°*>, respectively. 
The corresponding stellar masses at these peaks are 
log(M s /M©)= 10.30 ± 0.25 for blue galaxies and 
log(M s /M Q ) = 10.40 ± 0.25 for red galaxies. These 
masses are around 0.23 and 0.30 the characteristic 
stellar mass M* 1O 1O ' 93 M of the overall YMB09 
GSMF, respectively. The maximum difference be- 
tween the blue and red mean M s -Mh relations is 
attained at log(Mh/M Q )sa 11.9; at this mass, the / s 
value of the former is 1.7 times smaller than the / s 
of the latter. For larger masses this difference de- 
creases. 

At the low-mass end, roughly / s oc Mh (oc M s 0,5 ) 
and f s oc (oc M s °' 8 ) for the blue and red 

samples, respectively, while at the high-mass end, 
/„ cx M-° A (oc M-°- 7 ) and /„ oc M^ ' 6 (oc A/7 15 ), 
respectively. 

It is important to note that the differences be- 
tween blue and red M s -M^ relations at almost all 
masses are within the la uncertainty of our infer- 
ences. We conclude that the M s -M^ (/ S -Mh) re- 
lation does not depend significantly on galaxy color 
(type). If any, the mean / s ~Mh relation of red galax- 
ies is narrower and more peaked than the one of 
blue galaxies. In the mass range where the abun- 
dances of blue and red galaxies are closer (10.0 < 
log(M s /M©)< 10.7), the intrinsic scatter around the 
M s -Mh relation would slightly correlate with color 
in the sense that the redder (bluer) the galaxy, the 
larger (smaller) is its M s for a fixed Mh, with a max- 
imum average deviation from the mean due to color 
not larger than ~ 0.1 dex. For masses smaller than 
M s w 10 9 ' 7 Mq, the correlation of the scatter with 
color would invert. 

The (slight) differences between blue and red M s - 
Mh (/s-Mh) relations can be understood basically by 
the differences in the respective cumulative GSMFs 
and, at a minor level, by the differences of the cor- 
responding HMFs for each case. The sharp peak in 
the red / S -Mh relation is associated to the turn-over 
at M s - 10 10 - 5 M Q in the GSMF of red galaxies (see 
Fig. IT]). 

In order to estimate the influence of the correc- 
tions introduced to the HMF for blue (red) galaxies, 
we have redone the analysis by using the original 
distinct (S-T) HMF without any correction but re- 
normalised to obtain the same fraction of halos as the 
fraction implied by the GSMF of blue (red) galax- 
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TABLE 1 
FIT PARAMETERS 



1 -) J 

Parameter 


All 

All 


Blue 


Red 


log M ,/, 


11.97 


11.99 


11.87 


logM* 


10.40 


10.30 


10.40 


P 


0.34 


0.37 


0.18 


a 


1.45 


0.90 


1.50 


7 


0.90 


0.90 


0.90 


a (M s < M*) 


0.000 


0.125 


0.000 


a (M s > M*) 


0.095 


0.125 


0.093 



ies with respect to the total GSMF . The results are 
shown in Fig. |B]with solid curves of blue color (blue 
galaxies) and red color (red galaxies). For compar- 
ison, the corresponding relations plotted in Fig. [5] 
are reproduced here (dot-dashed blue and red lines, 
respectively). One sees that the corrections to the 
HMF we have introduced for associating halos to 
the blue and red galaxy sub-samples act in the di- 
rection of reducing the differences among them in 
the M s -Mh (/ S -Mh) relations, specially for larger 
masses. The group/cluster mass function correction 
to the HMF hosting central blue galaxies is the dom- 
inant one. The dashed blue and red curves show 
such when only this correction (and a small 

re-normalisation) is applied. 

3.2.1. Analytical fits to the stellar-halo mass 
relations 

From the comparison of the GSMF and HMF it 
is easy to deduce that high- and low-mass galaxies 
have significantly different M s -Mh scalings, a fact 
attributed to the different feedback/gas accretion 
mechanisms dominating in large and small systems 
(see e.g., Benson et al. 2003). The transition point 
between the low- and high-mass scalings defines a 
characteristic halo mass Mq^ and an associated stel- 
lar mass M* . Therefore, it was common to describe 
the M s -Mh relation as a double-power law with the 
turnover point at Mq^- However, BCW10 have ar- 
gued recently that a power-law at the high-mass side 
is conceptually a bad description for the M s -Mh re- 
lation and proposed a modification to it. Our results 
show indeed that a power-law is not enough to de- 
scribe the high-mass side of the M s -Mh relations. 

We have found that a good analytical descrip- 
tion to the overall, blue, and red mean M s -Mh re- 
lations inferred here can be obtained for the in- 
verse of the relations (Mh as a function of M s , as 
in BCW10) by proposing a power-law dependence 
for low masses and a sub-exponential law for high 
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Fig. 7. Analytical fits given by eq. ([5|) and Table 1 
compared to the mean M s -Mh relation obtained here 
for all central galaxies (black solid line) and the central 
blue (blue solid line) and red (red solid line) galaxy sub- 
samples. 

masses (see BCW10). The functional form that fits 
well the three Mh-M s relations is: 



M, 



0,h 



2"' 



Ms 
M* 



M« 
M* 



a/j 



1 T 



10 



a(M a /M*-l) 



(5) 

where (3 regulates the behaviour of the relation 
at masses M s < M*, a together with the sub- 
exponential term (a < 1) regulate the behaviour at 
masses M s > M*, and 7 regulates the transition of 
the relation around M* . In Table 1 are given the val- 
ues of all the parameters that best fit our results for 
the (central) overall, blue, and red M s -Mh relations. 
Note that a assumes two different values depending 
on whether the mass is smaller or larger than M*. 

Figure [7] shows the three mean M s -Mh relations 
obtained here and the functional form given in eq. 
([5]) with the corresponding parameters reported in 
Table 1. The functional form is an excellent fit to 
the overall and blue M s -Mh relations at all masses 
and to the red M s -Mh relation for masses larger than 
M h w 10 11 ' 3 M . 

3.3. The baryonic-halo mass relations for central 
blue and red galaxies 

The right upper and lower panels of Fig. [5] show 
the mean M^-M^ and /b-Mh relations, as in the left 
panels, for all central galaxies (solid line), central 
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blue (long-dashed blue line), and central red (short- 
dashed red line) galaxies. The blue and red GBMFs 
were calculated from the corresponding GSMFs and 
adding to M s the respective gas mass, M g (see §§2.1). 
The total GBMF is the sum of both of them. The 
error in Mb was calculated as the sum in quadratures 
of the errors in M s and M g . This error, together with 
the intrinsic scatter in M s (see §§2.2), both propa- 
gated to the M b -M h relation, account for an uncer- 
tainty (standard deviation) of ~ 0.23 dex in logMb 
at all masses (grey curves connected by vertical lines 
in Fig. |5j}. 

The baryonic mass fraction, fb, for blue galaxies 
is larger than the corresponding stellar one, / s , in 
particular at smaller halo masses. At Mh re 10 11 
M Q , /b is a factor 2.4 times higher than / s , while the 
peak of / b = 0.028±g;gif (at M h = 10 120 M ) is only 
1.3 times larger than the peak of / s (at Mh = 10 120 
Mq). For larger masses, the difference between / b 
and f s decreases, while for smaller masses, the lower 
is Mh, the larger is /b than /„. For red galaxies, 
f s and /b are very similar, some differences being 
observed only at the lowest masses. 

For masses larger (smaller) than Afh ~ 10 116 
M©, the differences between the M b -Mh (/b~Mh) 
relation of blue and red galaxies become smaller 
(larger) than in the case of stellar masses (see §§3.2 
and left panels of Fig. [5j. In general, the /b bell- 
shaped curve for red galaxies is more peaked and 
narrower than the one for blue galaxies. 

For blue galaxies, roughly f b oc M h ' 7 (Mb ' ) 
at the low-mass end, and fb oc Mh~ ' (Mb - " ) at 
the high-mass end. For red galaxies, roughly /b oc 
Mh 2 - 9 (M b - 8 ) at the low-mass end, and f h oc M h ~ - 6 
(Mb ) at the high-mass end. For halos of masses 
M h re 10 u o Mq and M h re 10 13 2 Mq, the baryon 
fraction for blue (red) galaxies decreases to values 
/b re 0.004 and 0.0085 (/ b re 0.0031 and 0.0071), 
respectively. Therefore, for all masses, /b << fu, 
where fu = ttb/Qni is the universal baryon mass 
fraction; for the cosmology used here, fu = 0.167. 

4. DISCUSSION 

4.1. Comparison with other works 

As discussed in §3.1 (see Fig. |4j, our inference 
of the local overall M s -Mh relation is in general in 
good agreement with several recent works that make 
use of the AMT (e.g., Baldry ct al. 2008; Guo ct al. 
2010; Moster et al. 2010; BCW10). The aim in this 
paper was to estimate the M s -Mh and Mb-Mh rela- 
tions for blue (late-type) and red (early-type) central 
galaxies separately. We have found that the differ- 
ences between the means of the obtained relations 



for blue and red galaxies are within the la uncer- 
tainty (see Fig. [5J. In more detail, the mean stellar 
and baryonic mass fractions (/ s and fb) as a func- 
tion of Mh for red galaxies are narrower and more 
peaked than those for blue galaxies in such a way 
that in a given mass range (11.5-13.0 and 11.5-12.5 
in log(Mh/M ) for the stellar and baryonic cases, re- 
spectively) the former are higher than the latter and 
outside these ranges, the trend is inverted, specially 
at the low-mass side. 

There are only a few previous attempts to infer 
the halo masses of central galaxies as a function of 
mass (luminosity) and galaxy type (Mandelbaum et 
al. 2006; More et al. 2011). These works use di- 
rect techniques (see Introduction), which are, how- 
ever, limited by low signal-to-noise ratios, specially 
for less massive systems in such a way that the halo 
mass determinations are reliable only for galaxies 
with M s £ 10 10 Mq. These techniques are galaxy- 
galaxy weak lcnsing and kinematics of satellite galax- 
ies around central galaxies. In order to overcome 
the issue of low signal-to-noise ratios in the current 
measures, large samples of galaxies are stacked to- 
gether in bins of similar properties (e.g., luminosity, 
M s , galaxy type) obtaining this way higher (statis- 
tically averaged) signals of the corresponding mea- 
sures (the tangential shear in the case of lensing and 
the weighted satellite velocity dispersion in the case 
of satellite kinematics). Besides, estimates of Mh 
with these sophisticated techniques are subject to 
several assumptions, among them, those related to 
the internal halo mass distribution. It is usual to 
assume the Navarro, Frenk & White (1997) density 
profile with the mean concentration for a given mass 
as measured in N-body cosmological simulations. 

It is not easy to achieve a fair comparison of the 
results obtained with the AM formalism and those 
with the direct methods. We have inferred the mean 
(and scatter) of logM s as a function of Mh, while the 
weak lensing and satellite kinematics techniques con- 
strain Mh as a function of M s (see e.g., More et al. 
2011); besides, the former calculates the mean of Mh 
(and its scatter) in a linear scale instead of a logarith- 
mic one. These different ways of defining the rela- 
tionship between stellar and halo masses, depending 
on the shapes and scatters of the corresponding rela- 
tions, diverge less or more among them. In BCW10 
(see their Fig. 10), it was shown that at low masses 
(log(M h /M Q < 12, log(M s /M Q ) < 10.5), averaging 
logAf s as a function of Mh or logMh as a function of 
M s give equivalent results for the AMT, but at high 
masses, where the M s -Mh relation becomes much 
shallower, this relation becomes steeper (higher stel- 
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lar mass at a fixed halo mass) for the latter case with 
respect to the former one. 

In Fig. the results from Mandelbaum et al. 
(2006) are reproduced, left panels for central late- 
type galaxies and right panels for central early-type 
galaxies (solid squares with error bars). The er- 
ror bars are 95 percent confidence intervals (statis- 
tical). Mandelbaum et al. (2006) have used the 
(de Vaucouleours/exponcntial) bulgc-to-total ratio, 
frac_deV, given in the SDSS PHOTO pipeline as 
a criterion for late- (frac_deV< 0.5) and early-type 
(frac_dV> 0.5) separation. This criterion of course is 
not the same as the color used in YMB09, but there 
is a correlation between both of them in such a way 
that a comparison between our results and those of 
Mandelbaum et al. is possible at a qualitative level. 
Note that we have diminished the halo masses of 
Mandelbaum et al. by ~ 15% on going from their 
to our definition of halo virial mass. In more recent 
works, Mandelbaum ct al. (2008) and Schulz et al. 
(2010) reported a new weak lensing analysis for the 
massive central early-type galaxies using the seventh 
SDSS data release (DR7) and a more sophisticated 
criteria for selecting the early-type lens population. 
Their results are plotted in the right panel of Fig. [8] 
with solid triangles and open squares, respectively. 

In the case of the satellite kinematics determina- 
tions of M h by More et al. (2011), the same SDSS 
sample and similar recipes as in YMB09 for cal- 
culating Ms, classifying galaxies into blue and red, 
and finding central and satellites galaxies were used. 
More ct al. (2011) applied their analysis to con- 
strain the mean logMh as a function of M s , but also 
present the constraints of their model for the mean 
of logM s as a function of Mh. Their results for the 
latter case, kindly made available to us in electronic 
form by Dr. S. More, are reproduced in Fig. [8] as 
the shaded (orange) regions which represent the 68% 
confidence intervals. On going from their to our def- 
initions of halo mass and IMF, their Mh and M s 
were diminished by w 15% and « 25%, respectively. 
The dotted horizontal lines in each panel show the 
approximate range in M s , where the determinations 
are reliable according to More et al. (2011; see their 
Fig. 11). 

In More et al (2011) are also reported results for 
the average Mh as a function of M s split in central 
blue and red galaxies corresponding to the galaxy 
group analysis by Yang et al. (2007). The solid 
(cyan) curves in Fig. [8] reproduce these results. 

Finally, the standard ±ler deviation intervals 
that we have obtained from the AMT are reproduced 
in Fig. [8] for central blue and red galaxies (solid blue 



and red curves connected by vertical lines, respec- 
tively). Note that in the determinations with direct 
methods, the systematic uncertainty in M ai which is 
the main source of error in the AMT, was not taken 
into account. 

Our inference for early-type (red) galaxies is con- 
sistent (within the uncertainties, errors, and different 
ways of presenting the constraints) with the weak 
lensing results of Mandelbaum ct al. (2006) and 
Schulz ct al. (2010), and with the galaxy group anal- 
ysis of Yang et al. (2007) as reported in More et al. 
(2011), for all the masses reported in each one of 
these papers. With respect to the satellite kinemat- 
ics analysis by More et al. (2011), their mean halo 
masses for M s ~ 5 x 10 9 — 10 11 M© (for smaller 
masses their uncertainties are very large), are larger 
than ours (and those of Mandelbaum et al. 2006) by 
factors around 2. For larger masses, all determina- 
tions agree roughly with our results. In fact, there is 
some indication that satellite kinematics yields halo 
masses around low mass central galaxies that are 
systematically larger than most other methods, spe- 
cially for red central galaxies (Skibba et al. 2011; 
but see More et al. 2011 for a discussion). 

For late-type (blue) galaxies, our results are in 
reasonable agreement with those of Mandelbaum et 
al. (2006) for masses M s < 10 10 8 M Q . At higher 
masses, their results imply halo masses for a given 
Ms smaller than ours, with the difference increas- 
ing to higher stellar mass. The discrepancy would 
be weaker if taking into account that the mean M s - 
Mh relation in our case becomes steeper when cal- 
culating Mh as a function of M s . On the other 
hand, it must be said that the number statistics be- 
comes poor for massive late-type galaxies, resulting 
in a stacked weak lensing analysis with large error 
bars. For example, in the two most massive bins 
in the Mandelbaum et al. (2006) sample (the two 
uppermost points in Fig. [8|), only 5 and 11 percent 
of the galaxies are classified as late types. Future 
weak lensing works should confirm whether high- 
mass late-type galaxies have or not such relatively 
small halos as found in Mandelbaum et al. (2006). 
Regarding the comparison with the satellite kine- 
matics inferences of More et al. (2011), the agree- 
ment is reasonable at least up to M s w 10 11 M Q , 
though the relation inferred by these authors is less 
curved than ours. For larger masses, these authors 
caution that their results become very uncertain, as 
in the weak lensing case, because of poor statistics 
of massive blue galaxies. The galaxy groups (Yang 
et al. 2007) inference, in the mass range allowed 
by this technique, gives halo masses slightly smaller 
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Fig. 8. Comparison with other observational inferences. Left panel: M s — Mh relation for blue (late-type) galaxies. 
The blue curves connected by vertical lines encompass the ±1<t interval inferred here. We also reproduce the inferences 
using galaxy-galaxy weak lensing by Mandelbaum et al. (2006, black squares), galaxy groups (Yang et al 2007, cyan 
solid line), and satellite kinematics (More et al 2011, orange vertical lines). Estimates for the Milky Way are plotted 
(open circle with error bar). Right panel: M s — Mh relation for red (early-type) galaxies. The red curves connected 
by vertical lines encompass the ±lcr interval inferred here. Other determinations as in the left panel but for early- type 
galaxies are plotted. More recent inferences with the weak lensing technique by Mandelbaum et al. (2008, filled violet 
triangles) and by Schulz et al. (2010, open green squares) are also plotted. 



than the means of our inference for a given M s . 

In general, most techniques for inferring the rela- 
tionship between stellar and halo masses of galaxies 
agree among them within factors up to 2-3 in Mh 
(BCW10; More et al. 2011; Dutton et al. 2010). 
This seems to bo also the case for samples parti- 
tioned into late- and early-type galaxies, as shown 
here. However, beyond the detailed comparison be- 
tween our results and those obtained with direct 
techniques, it seems that there is a systematic qual- 
itative difference: in our case, at a given halo mass 
(for 10 115 Mq < M h < 10 130 M Q ), blue centrals, on 
average, have lower stellar masses than red centrals, 
while in the case of determinations with direct tech- 
niques, the opposite applies at least for masses larger 
than M h ~ 10 12 M (Mandelbaum et al. 2006; More 
et al. 2011; see also Figs. Eland©. 

A partial source of bias contributing to this dif- 
ference could be that in the weak lensing and satellite 
kinematics techniques the same concentration for ha- 
los hosting late- and early-type galaxies is assumed. 
If halos of late- (early-)type galaxies are less (more) 
concentrated than the corresponding average, then 
for the same measure (shear or satellite velocity dis- 
persion), the halo masses are expected to be higher 
(lower) than the obtained ones. Therefore, the dif- 



ferences found (Mandelbaum et al. 2006 and More 
et al. 2011) in the mass halos of late- and early-type 
galaxies of a given M s would decrease or even invert 
their sense. 

While it is difficult to make any robust state- 
ment about possible systematics in one or another 
technique regarding late and early types, we ask our- 
selves what should be modified in our assumptions in 
order to invert the behaviour of the M s -Mh relations 
with galaxy type (color) obtained here. We have 
shown in Fig. |6] that our corrections to the HMF 
had the effect of making closer the M s -Mh relations 
of blue and red galaxies at large masses. One possi- 
bility in order not only close more the relations but 
invert them is to make even steeper (shallower) the 
HMF corresponding to blue (red) galaxies, mainly 
at the high-mass end (see Fig. [3j lower panel). This 
would imply, for instance, a higher correction to the 
HMF due to groups than that made by us. The 
group/cluster mass function used by us (Hcinamaki 
et al. 2003) is one of the most general ones found in 
the literature; it includes all kinds of groups/clusters 
with 3 or more members and SN/N > 80. The au- 
thors note that their sample is complete down to a 
dynamical mass roughly equivalent to Mh = 5 x 10 13 
M Q . It could be that the abundance of groups of 
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lower masses is larger than that given in Heinamaki 
et al. (2003), though it is difficult to accept that 
blue galaxies are completely absent in the centers of 
small and loose groups of a few (> 2) members. 

Last but not least, in Fig. [8] we include observa- 
tional estimates for our Galaxy (open circle). The 
uncertainties in the estimates of Mh for the Milky 
Way are still large but better than most of the deter- 
minations for other individual galaxies. For recent 
reviews on different results see Guo et al. (2010) 
and Dutton et al. (2010). In Fig. [8] we plot a recent 
estimate of Mh based on observations of 16 high ve- 
locity stars (Smith et al 2007). These authors find 
M h = 1.42±o.54 x 10 12 M Q , which is in good agree- 
ment with several previous works (e.g., Wilkinson 
& Evans 1999; Sakamoto et al. 2003; Li & White 
2008), though results from Xue et al. (2008) suggest 
lower values (but see a recent revision by Przybilla 
et al. 2010). For its M h , the M s of Milky Way seems 
to be in the high extremum of blue galaxies, close to 
values typical of red galaxies. It should be said that 
it is an open question whether the Milky Way is an 
average galaxy or not. In the stellar Tully-Fisher and 
radius-M s relations (e.g., Avila- Reese et al. 2008), 
the Milky Way is shifted from the average to the 
high-velocity and low-radius sides, respectively. 

4.2. Interpretations and consistency of the results 

Although our main result is that the differences 
between the M s -Mh and Mb-Mh relations for cen- 
tral blue and red galaxies are marginal (within the 
uncertainties of our determinations), we will explore 
whether such differences are expected or not. For 
this it is important to approach the problem from 
an evolutionary point of view. 

In Firmani & Avila- Reese (2010, hereafter FA10), 
the determinations of the M s -Mh relation for all 
galaxies at different redshifts, out to z = 4 (BCW10), 
and the average ACDM individual halo mass ag- 
gregation histories (MAHs) were used to determine 
the individual average M s growth of galaxies in gen- 
eral as a function of mass (called in that paper 
as Galaxian Hybrid Evolutionary Tracks, GHETs). 
It was found that the more massive the galaxies, 
the earlier transit from their active (star-forming, 
blue) regime of M s growth to a passive (red) phase 
(population 'downsizing'), while their correspond- 
ing halos continue growing, more efficiently at later 
epochs as more massive they are ('upsizing'). The 
inferred trend for the transition stellar mass is 
log(M tran /M©)« 10.30 + 0.55z. Therefore, galax- 
ies of mass M s w 10 10 3 M© are on average becom- 
ing passive (red) today. For M s ^ M tran , the larger 



the mass, the redder will be the galaxy on average. 
The opposite applies for M s ^ M tran , the smaller 
the mass, the bluer will be the galaxy. Interesting 
enough, M s w 10 10 3 M Q is roughly the mass where 
the overall YMB09 blue and red GSMFs cross: for 
masses larger than this crossing mass, M cross , redder 
galaxies become more and more abundant than bluer 
ones and the inverse happens at smaller masses (see 
Fig. P). 

Galaxies that are transiting from active to pas- 
sive at z ~ (those around M tran sa 10 10 3 M ) have 
probably been subject recently to a process that in- 
duced an efficient transformation of the available gas 
into stars in such a way that their stellar popula- 
tions started to redden passively. Hence, for a given 
Mh, they are expected to have a higher M s (or / s ) 
than those galaxies of similar mass that did not suf- 
fer (yet?) the mentioned above process (bluer ones). 
The relatively small difference in / s for blue and red 
galaxies we have found here (whose maximum is at- 
tained namely around M tran ~ M cross , Fig. [SJ would 
imply that the scatter around M tran is moderate. 

Galaxies more massive than M tran (or M cross ), 
according to the evolutionary analysis by FA10, had 
the process of efficient gas consumption into stars 
(and the further cessation of M s growth) earlier on 
average as more massive is the galaxy, while their 
halos continue growing. Therefore, one expects that 
the more massive the galaxy, the redder and the 
lower its stellar (and baryonic) mass fraction / s will 
be on average. The few blue massive galaxies may 
have slightly smaller stellar masses (lower / s ) than 
the corresponding red ones because they should have 
transformed gas into stars less efficiently in the past. 
Therefore, by including gas, i.e. when passing to /b 
the difference between blue and red massive galax- 
ies at large masses should become negligible. This is 
what indeed happens (see Fig. [5]). 

Galaxies less massive than M tran (or M cross ) at 
z ~ 0, according to FA10, are in general more ac- 
tively assembling their stellar masses as smaller they 
are ('downsizing in specific SFR), while their dark 
halo mass growth is already very slow. This implies 
the existence of relatively larger reservoirs of cold gas 
in the galaxies as smaller they are (gas not related 
to the halo-driven infall) because the SF has been 
delayed in the disk and/or cold gas is being lately 
(re)accreted into the galaxy. However, if for some 
reason the gas reservoir in these galaxies is lost, then 
the galaxy will redden and its baryonic and stellar 
mass fractions will be smaller than of the galaxies 
that were able to keep their gas reservoir (the ma- 
jority), in agreement with our inferences here (Fig. 
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0). 

5. SUMMARY AND CONCLUSIONS 

By means of the AM technique and using the 
central blue and red GSMFs, constructed from the 
local SDSS sample by YMB09, we have inferred the 
local M s -Mh (or / S -Mh) relations for central galaxies 
and for the sub-samples of blue and red galaxies. To 
derive the relations for the sample of blue galaxies, 
(i) the mass function of observed groups/clusters of 
galaxies is subtracted from the distinct (S-T) HMF 
(blue, late-type galaxies are not observed in the cen- 
ters of groups and clusters), and (ii) halos that suf- 
fered a major merger since z = 0.8 are excluded. For 
red galaxies, the HMF is assumed to be the comple- 
ment of the "blue" one, with respect to the overall 
(distinct) HMF . We consider as sources of uncer- 
tainty in our analysis only the systematical error in 
assigning stellar masses to galaxies (0.25 dex) and 
the intrinsic statistical scatter in stellar mass at a 
fixed halo mass (0.16 dex). By using the observa- 
tional Mg-M s relation and its scatter, we transited 
from M s to M b (=M S + M g ) in the GSMF and es- 
timated the overall blue and red GBMFs, which are 
used to obtain the corresponding baryonic M^-M^ 
(or /b~Mh) relations using the AM technique. 

The M s -Mh relation obtained here agrees rather 
well with previous studies (see Fig. H]). The small 
differences found in this work can be explained 
mainly in terms of the different GSMFs used in 
each study, and to a less extent by variations in the 
methodology. The la uncertainty in the obtained 
M s -M h relation is « 0.25 dex in logM s . The Ms- 
Mh relation of central galaxies lies below (lower M s 
for a given Mh) the overall one by a factor ~ 1.6 at 
M h = 10 11 M Q and less than 5% for M h > 10 13 M . 

Our main result refers to the calculation of the 
central M s -Mh and Mb-Mh relations for the two 
broad populations into which the galaxy sample can 
be divided: blue (late-type) and red (early-type) 
galaxies. We highlight the following results from our 
analysis: 

• At Mh ~ 10 113 M Q the mean stellar mass frac- 
tion f s of blue galaxies is lower than the one of red 
galaxies, the maximum difference being attained at 
M h w lO 117 M ; at this mass, the / s of red galaxies 
is 1.7 times the one of blue galaxies (see Fig. [5]). At 
larger masses, the difference decreases until it disap- 
pears. At Mh ^ 10 113 M Q the trend is reversed as 
blue galaxies tend to have higher values of / s than 
red ones. In the case of the baryonic mass fractions, 
/b , the same trends of the stellar relations remain but 
at Mh ^ 10 113 M Q the difference in /b between blue 



and red galaxies is small, while for smaller masses, 
the difference increases. 

• The M s -M h and M h -M h (or f B -M h and f b - 
Mh) relations of central blue and red sub-samples do 
not differ significantly from the respective relations 
of the overall central sample, and these differences 
are within the la uncertainty of the inferences (Fig. 
[5]). For blue (red) galaxies, the maximum value of / s 
is 0.021+° °^ (0.034+°-°i5) and is attained in halos 
of mass M h = 10 1198 M Q (M h = 10 1187 M Q ); the 
corresponding stellar mass is M s = io 10 - 30±0 - 25 M 
(M B = 10 10 - 40±0 - 25 Mq), which is around 0.23 (0.30) 
times M*, the Schechter fit characteristic mass of 
the overall GSMF of YMB09. For smaller and larger 
masses, / s significantly decreases. 

• We have compared our results with the few ob- 
servational inferences of the M s -Mh relation for blue 
(late-type) and red (early-type) galaxies that exist in 
the literature. Although these studies estimate halo 
masses using direct techniques (weak lensing and 
galaxy satellite kinematics), they are still limited by 
the stacking approach they need to apply (due to the 
low signal-to-noise ratio of individual galaxies) and 
by the large uncertainties owing to the unknown sys- 
tematics. The overall differences among the different 
studies (including ours) amount up to factors 2-4 at 
a given mass (these factors being much smaller at 
other masses) for most methods (Fig. [8]). For blue 
galaxies, all methods agree reasonably well for low 
masses (Mh £3 x 10 12 M ), but at higher masses, 
our inference implies larger halos for a given M s than 
the results from direct techniques. For red galax- 
ies, at high masses (Mh ^ 3 x 10 12 M Q ), all methods 
agree reasonably well, but at lower masses, the satel- 
lite kinematics technique produces halo masses, for a 
given M s , larger than those obtained by other meth- 
ods. 

• According to our results, for Mh 10 11 3 Mq, 
the intrinsic scatter of the M s -Mh relation should 
slightly anti-correlate with galaxy color (for a fixed 
Mh, the bluer the galaxy, the higher its M s ), while 
for more massive systems, the correlation should be 
direct (for a fixed Mh, the redder the galaxy, the 
higher its M s ). For massive blue galaxies in order 
to had have higher higher / s values than the red 
ones as the results from direct techniques suggest, 
the HMF halos hosting blue (red) galaxies should be 
even steeper (shallower) than what we have proposed 
here; this seems too exaggerated. 

• The maximum baryon mass fraction of blue 
and red galaxies arc /b = 0.028^q'q^ 8 and /b = 
0.034^q q 2 4, respectively, much smaller than fu = 
0.167 in both cases, and these maxima are attained 
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at Mh « 10 12 M Q . At large masses /b decreases ap- 
proximately as fb oc Mh _0 ' 5 (Mb -0 ' 8 ) for blue galax- 
ies and as /b oc A/h"' 6 (^b _1 ' 5 ) for red galaxies, 
in such a way that from Mh « 5 x 10 12 M©, blue 
galaxies have on average slightly larger values of 
/b than red ones. At low masses, the /b of red 
galaxies strongly decreases as the mass is smaller 
/ b oc M h 2 - 9 (M b - 8 ), while for blue galaxies, due to 
the increasing gas fractions as smaller is the mass, /b 
decreases slower than / s , as /b oc Mh ' 7 (oc Mb ' ). 

The AM technique has been revealed as a rel- 
atively simple but powerful method for connecting 
empirically galaxies to dark halos. Here we extended 
this technique towards inferences for the blue and red 
galaxy sub-populations separately. By introducing 
a minimum of assumptions -otherwise the method 
becomes close to a semi-analytical model™ we have 
found that the stellar and baryon mass-halo mass 
relations of blue and red galaxies do not differ sig- 
nificantly among them and from the overall ones. 
The maximum differences are around the peak of 
these relations, Al^ w 10 12 M , and are consistent 
qualitatively with the inference that the galaxies in 
these halos are transiting from an active to a quies- 
cent regime of M s growth (FA10). Those that tran- 
sited recently is because they had an efficient process 
of gas consumption into stars and further cessation 
of Ajf s growth; therefore, should be redder and with 
higher / s values than those that still did not transit. 
For larger and lower masses than Mh ~ 10 12 M , 
the differences decrease and even invert, something 
that is also consistent with the inferences by FA10, 
based on the semi-empirical determinations of the 
evolution of the overall M s -Mh relation. 
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